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ABSTRACT 

We analyze 40 cosmological re-simulations of individual massive galaxies with present-day stellar 
masses of M* > 6.3 x lO M© in order to investigate the physical origin of the observed strong 
increase in galaxy sizes and the decrease of the stellar velocity dispersions since redshift z k, 2. At 
present 25 out of 40 galaxies are quiescent with structural parameters (sizes and velocity dispersions) 
in agreement with local early type galaxies. At z=2 all simulated galaxies with M* > 10 n M© (11 
out of 40) at z=2 are compact with projected half-mass radii of « 0.77 (±0.24) kpc and line-of-sight 
velocity dispersions within the projected half-mass radius of w 262 (±28) kms -1 (3 out of 11 arc 
already quiescent). Similar to observed compact early- type galaxies at high redshift the simulated 
galaxies are clearly offset from the local mass-size and mass-velocity dispersion relations. Towards 
redshift zero the sizes increase by a factor of ~ 5 — 6, following R\/2 tx (1 + z) a with a = —1.44 for 
quiescent galaxies (a = —1.12 for all galaxies). The velocity dispersions drop by about one-third since 
z rs 2 , following ai/ 2 oc (1 + z)" with (3 — 0.44 for the quiescent galaxies (/? = 0.37 for all galaxies). 
The simulated size and dispersion evolution is in good agreement with observations and results from 
the subsequent accretion and merging of stellar systems at z < 2 which is a natural consequence of 
the hierarchical structure formation. A significant number of the simulated massive galaxies (7 out of 
40) experience no merger more massive than 1:4 (usually considered as major mergers). On average, 
the dominant accretion mode is stellar minor mergers with a mass- weighted mass-ratio of 1:5. We 
therefore conclude that the evolution of massive early-type galaxies since z ~ 2 and their present-day 
properties are predominantly determined by frequent 'minor' mergers of moderate mass-ratios and 
not by major mergers alone. 

Subject headings: galaxies: elliptical - galaxies: interaction- galaxies: dynamics - galaxies: evolution 
- methods: numerical 



1. INTRODUCTION 

There is growing observational evidence for the exis- 
tence of a population of massive galaxies (w 10 11 Mq) 
with small sizes (~ 1 kpc) and low star for- 
mation rates at redshift z > 2. These galax- 
ies are smaller by a factor of three to five com- 
pared to present-day e llipticals at similar masses (e.g. 
IHvde fe Bernard! I2009f) and their effective stellar den- 
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2009; Bczanson ct al. 2009; Damianov ct al. 2009, 2011). 



Deep observations down to low surface brightness limits 
(H w 28 mag arcsec -2 ) show no e vidence for faint, pre- 
viously missed , stell ar envelopes (|Carrasco et al.1 120101 : 
ISzomoru et al.1 120101 ) and measurements of higher ve- 
locity dispersions seem to in dependently confirm pre- 
vious high mass estim ates (Ivan der Wei et al~l 1 20051 
2009HCenarro fc Truiillol 120091; Ivan Dokkum et al.l|2009l: 
Cappellari et al.ll2009t iMartinez-Manso et al.ll2011f ). 



Quiescent (red & dead) galaxies make up about half of 



the general high redshift (z k. 2) population of massive 
galaxi es and most of them (90%) are found to b e com- 
pact dKriek et all [2001 Ivan Dokkum et al.l [2001 12001 
I Williams et al.ll2009f ). In the local Univer se, galaxies of 
simil ar mass and size are ext remely rare (iTruiillo et al.l 
120091) or do not exist at all (jTavlor et all I2010D . This 
indicates that present-day early-type galaxies were not 
fully assembled at z » 2 and underwent significant struc- 
tural evolution until the present day. Observations of the 
growth of massive gal axies since z ~ 2 selecte d at con- 
stant number density ([van Dokkum et aLll2010f) indicate 
that they grow inside-out. A quiescent - without signifi- 
cant in-situ formation of new stars - build-up of extended 
stellar envelopes can originate from minor mergers and 
was predicted from cosmologi cal simulations (jNaab et al.l 
[2001 [20091 IQser et al-lllHOl ) and recently, for the first 
time such minor merge rs at high redshift might have 
been directly observ ed (jvan Dok kum fc Brammerl [2010; 
iCarrasco et^l2010h . 

A simple picture of high redshift monolithic formation 
or, similarly, a binary merger of massive very gas-rich 
disks at z > 2 - which has been suggested as a rea- 
sonable formation mechanism for compact high-redshift 
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galaxies (iWuvts et all 120101 iRicciardelli et al.l 120101: 
iBournaud et al.l 1201 If) - followed by passive evolution 
can be ruled out (iKriek et al.l 120081 : Ivan Dokkum et al] 
120081: iTruiillo et aLf TiToOQ'l unless the increase in size of 
ellipticals can be explained by secular processes such as 
adiabatic expans ion driven by stellar mass loss and/or 
strong feedback (jFan et al.l 120081: iDamianov et al.l 120091: 
iFan et all I2010D . T his process seems to be disfavored 
by ob servations (e.g. iHopkins et al.|[2010bl : ITruiillo et al.l 
1201 If ) and the absence of a significant young stellar pop- 
ulation would indicate that such hypothesized secular 
processes would need to occur without significant star 
formation. 

Based on high-resol ution cosmo l ogical simul ations 
of individual ga laxies, iNaab et all (|2007l 12009ft and 
IQser et all (|2010D provide an explanation for the size 
growth and the decrease in velocity dispersion, which 
is consistent with the cosmological hierarchical buildup 
of galaxies. The compact cores of massive galaxies 
form during an early rapid phase of dissipational in- 
situ star formation at 6 > z > 2 fed by cold flows 
(|Keres et al.l l200l IDekel et al.l 12001 lOser et al.l [20101) 
and/or gas ric h mergers leading to large stellar sur- 
face densities (jWeinzirl et al.l 120111 ). At the end of 
this phase the observed as well as simulated galaxies 
are more flattened and disk-like than their low red- 
shift counterparts (Ivan Dokkum et al.l [200 8: Naa b et ail 
20091 Ivan der Wei et al.ll201lD . Thevare already massive 
(« IO^Mq) but have small sizes of « 1 kpc and velocity 
dispersions of « 240 kms" 1 (see also lJoung et al.H2009l) . 
in general agreement with observations. The subsequent 
evolution is dominated by the addition of star s that have 
forme d ex-situ, i.e. outside the galaxy itself (Oscr et al. 
120101) . T hese accreted stars ty pically settle at larger radii 
(see also lCoccato et af1l2011l) . The early domination of 
in-situ star formation and the subsequent growth by stel- 
lar mergers is in ag reement with pred i ctions from semi- 
analyt i cal models (Kauffm annl Il996t IKhochfar fe Silkl 



2006 at iDe Lucia et all 120061; iDe Lu cia & Bl aizotl 1200 
Guo fc White! 120081: IShankar et al.ll2010allbD and the as 



sembly scenario discussed in IKormendv et all (|2009T) ■ 

In the absence of gas, stellar (i.e. collisionless) ac- 
cretion and 'dry' merging in general is an energy con- 
serving process in the sense that none of the gravita- 
tional and binding energy in the accreted systems can 
be radiated away during the merging event. Therefore, 
while the galaxies grow in mass, they must significantly 
increase their sizes and, eventually, decrease their ve- 
locity dispersions during this phase, in particular if the 
stars are accreted in minor mergers. In massive galax- 
ies that are embedded in a hot gaseous halos some frac- 
tion of the gravitationa l energy can be radiated away 
Johansson et al.ll2009bfl. iCole et all (|2000T ). INaab et all 



2009) and IBezanson et al.l (|2009D presented the simple 



virial arguments for why minor mergers lead to a stronger 
size increase and a decrease in velocity disp ersion than 
the m ore commonly studied major mergers (INaab et al.l 
200C 



In lOser et al.l (|2010f ) we investigated this two phase 
scenario in more detail with a larger sample of re- 
simulations and found a connection bet ween galaxy mass, 
size, a n d the assembly history (see also IKhochfar fc Silkl 
l2006at IDe Lucia fc Blaizotl 120071: IGuo fc White! 120081: 



IHopkins et aTll2009bl: IFeldmann et al.ll2010l ). More mas- 
sive present-day systems contain a larger fraction of ac- 
creted stars (up to 80 per cent) which, over time, build 
an outer envelope and increase the size of the systems 
(INaab et al.l[2009t IHopkins et al.ll2009aT : lOser et aL|[20Tol : 
IFeldmann et al.1 120101 ) . This scenario receives support 
from recent observational findings that massive galax- 
ies have increased their mass at radii r > 5 kpc by a 
factor of tts 4 since z=2 with the mass at smaller rad ii 
being essentially unchanged (jvan Dokkum et aII[20T0| ). 
It is als o the favored mode l to explain observed kine- 
matics (lArnold et al.l 1201 lh and metallicity gradients 
(|Forbes et al.ll201lD of globular cluster populations in 



nearby elliptical galaxies. 

In this paper we analyz e a subset of massive galaxies 
from the lOser et all ()2010f ) simulations with a particular 
focus on the evolution of sizes and velocity dispersions 
since z « 2. The paper is organized as follows: In section 
[2] we briefly review the simulations. The results on the 
evolution of size and velocity dispersion are presented in 
sections [3] and [4] The stellar merger histories of the res- 
imulated galaxies are reviewed in section [5j We conclude 
and discuss our results in section |6l 

2. HIGH RESOLUTION SIMULATIONS OF INDIVIDUAL 
GALAXY HALOS 

The results presented in this paper are drawn from 
40 'zoom-in' hydrodynamic simulatio ns of individual ha - 
los which are presented in detail in I Oser et al.l (2010). 
The halos are picked from a dark matter only simula- 
tion using a fla t cosmology with p arameters obtained 
from WMAP3 (ISpergel et al.l 12001 : h = 0.72, ft h = 
0.044, n dm = 0.216, Q\ = 0.74, a 8 = 0.77 and an ini- 
tial slope of the power spectrum of n s = 0.95. From 
rcdshift zero we trace back in time all particles close 
to the halos of interest at any given snapshot. Those 
particles are then replaced with high-resolution gas and 
dark matter particles. The original dark matter par- 
ticles are merged (depending on their distance to the 
re-simulated halo) to reduce the particle count and the 
simulation time. The new high mass resolution initial 
conditions are evolved from redshift z=43 to the present 
day using a modi fied version of the parallel TreeSPH 
code GADGET-2 (ISpringell pOOl including star f prma- 
tion, supernovae feedback f Springcl fc Hernquistl 12003?) 
and cooling for a primordial composition of hydrogen and 
helium. Additionally, the simulations include a redshift- 
dep endent UV background radiation field with a modi- 
fied lHaardt fc Madaulll996l spectrum. 

The simulated halo masses cover the range 7 x 
1O 11 M /i~ 1 < M vir < 2.7 x lO^M©/*" 1 and the cen- 
tral galaxy masses are between 4.5 x 1Q W M^h^ 1 < 



0. The masses for 



4.2 x 1O 6 M /i" 1 



M* < 3.6 x IO^Mq/i- 1 at z 
the gas and star particles are m*. gas 
(we spawn one star particle per gas particle), whereas 
the dark matter particles have a mass of m^m — 
2.5 x 10 7 Mq/i -1 . The comoving gravitational soften- 
ing lengths used are e gas ,star = 400pch~ 1 for the gas 
and star particles and eh a io = 890pch~ 1 for the dark 
matter. At z w 2 the corresponding physical soften- 
ing length for stars is e gaS:Sta r = ISSpch" 1 . The inte- 
gration accuracy parameter wa s set to 0.005 t o guaran- 
tee accurate time integration (|Springelll2005T) . In the 
following we present the results for 40 galaxies with 
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masses larger than w 6.3 x 1O 1O M0 for direct com- 
parison with observations. These galaxies are well re- 
solved with w 1.5 x 10 5 — 3 x 10 6 particles within the 
virial radius (i? v ir = -R20O1 the radius where the spheri- 
cal overdensity drops below 200 times the critical den- 
sity of the universe at a given redshift). Using the 
above simulation parameters for zoom simulations have 
been shown to result in galaxies with reasonable present- 
day properties (INaab et ai1l2007HJohansson et al.f 2009b; 
INaab et all 120091 : IQser et all 120101 ). However, the frac- 
tion of available baryons converted into stars, /*, for 
galaxies in this mass range is typically 2 times higher 
than estimates from models that are constructed by 
matching observed luminosity functions to simulated 
halo mass functions (iGuo et al.ll2010l iMoster et aT]|2010t 
iBehroozi et alJl2010L and references therein). 

3. REDSHIFT EVOLUTION OF SIZES 

In the left panel of Fig. Q] we show the projected half- 
mass radii of the simulated galaxies as a function of their 
stellar mass at z=0 (blue circles), as well as the sizes 
and masses of their most massive progenitors at z=l 
(red triangles), z=2 (green triangles), and z=4 (purple 
stars). We determine the center of the g alaxies using 
the sh rinking sphere technique described in iPower et al.l 
(2003) , starting at redshift with a sphere that contains 
all the stellar particles. In all the previous snapshots 
the center of the most massive progenitor is found with 
the same technique starting with a sphere of a radius 
of at least 25kpc that encompasses the 50 innermost 
particles identified in the last processed snapshot. The 
sizes indicated here are the mean values of the half-mass 
radii of all stars within 0.1 x i? vir = i?Gai (considered the 
central galaxy) projected along the three principal axes 
of the main stellar body. We always show the median 
of the sizes of the galaxies that we compute from the 
snapshot at the target redshift as well as the two snap- 
shots before and after this one to avoid outliers caused 
interacting systems. We separate the sample into qui- 
escent galaxies (solid symb ols) with specific star forma- 
tion rates sSFR < 0.3/t H (|Franx et a ll 1200 8). where i H 
is the age of the Universe at each redshift. Star form- 
ing galaxies are indicated by open symbols and have 
sSFR > 0.3/t H . The dashed black line shows the z=0 
line ar fit to the SLACS sample of local early-type galax- 
ies (jNipoti et al.ll2009T ) with the observed scatter given 
by the dotted lines, which is in good agreement with the 
simulated galaxy sizes. Other published local mass-size 
relations have slightly d i fferent slopes and offsets (see e.g. 
Hyde fc Ber"nardil [200l IGuo et all [20091: iWilliams et all 
2010) which does, however, not affect our general conclu- 
sions. At z=4 all progenitor galaxies are very compact 
with similar sizes (~ 300 pc) independent of their mass. 
During this phase the formation of the proto-galaxies is 
dominated b y gas dissipation an d in-situ star formation 
(|Naab et all 120091: Uoung et all 120091: [OleTet^ [2010l ). 
By z=2 a clearly visible mass-size relation has already 
developed. At this epoch the most massive galaxies of 
our sample have experienced the most rapid size growth 
with half- mass radii up to rs 1.3 kpc for galaxies with 
10 in stellar mass, in good agr eement even with 

the most recent observations (e.g. iRvan et all 120101 : 
Ivan der Wei et all 120111: iCassata et al.l 120111 ). Towards 
z=0 the simulated galaxies continue to grow in size as 



well as mass. The descendants of galaxies that are al- 
ready massive (> 6.3 x 1O 1O M0) at z=2 (green symbols 
to the right of the vertical dashed line) increase their 
mass by a factor of 3.5 (77 ± 10% of the accumulated 
mass is due to stellar accretion) and their projected half 
mass radii grow by a factor of 6.5. On average, all sim- 
ulated galaxies more massive than 6.3 x 10 10 M Q at any 
given redshi f t grow by a factor of 2.1 in mass (see e.g. 
iFaber et~aT1 ([2007D ) and a factor of 4.6 in radius since 
z=2. This already indicates that the size growth cannot 
be the result of equal-mass dry mergers, since the ensuing 
size growth should be, at most, direct ly proportional to 
the mass increase (jNipoti et all 120031 120091: INaab et all 
2009). Overall, the size growth is differential, i.e. the 
most massive galaxies show the strongest size increase 
and the tilt of the mass-size relation increases towards 
lower redshifts. Fitting all galaxies with R1/2 cx we 
find 5 w 0.46 ±0.056 at z=2 and S » 0.69 ±0.064 at z=0. 
For quiescent galaxies (with worse statistics) we find no 
trend for differential size growth with S « 0.67 ±0.069 at 
z=2 and 5 ps 0.65 ± 0.090 at z=0. This i s in qualitative 
agreem ent with recent observations by IWilliams et all 
(2010) who do not find observational indications for a 
differential size growth of quiescent galaxies. 

Observed sizes of massive galaxies are found to evolve 
as (1 + z) a . Depending on the selection criteria (specific 
star formation rate, concentration etc.) and observed 
redshift range the observed val ues of a for massive el- 
lipticals ra nge from a = -0-75 (iNewman et~ani2010h to 
a = -1.62 (jDamianov et all 120111 ). iFranx et alT (2008 > l 
find a —0.8 for all galaxies in this mass range (M* > 
10 11 Mq) whereas for quiescent massive galaxies the ob- 
served size ev o lution is faster with —1.09 < a < —1.22. 
ICassata et all (|2011[ ) obtain values fr om a = -0.87 to a 
= -1.4 2 depending on stellar mass. Ivan der Wei et all 
(2008) find a value of a w -0.98, when they include re- 
sults from previous surveys this changes to a w -1.20. A 
similar trend ( a « —1.11) is found for UV-bright galaxies 
(jMosleh et al.ll20Tll) . 

In the right panel of Fig. [1] we show the size evolu- 
tion of galaxies m ore massive than 6.3 x 10 10 M© (see 
IFranx et al.l l200l since z=2. At z=4 all progenitor 
galaxies drop below the threshold mass, but are still re- 
solved by ~ 10 4 particles. We added the results from 
INaab et all (|2009| ). which were obtained with the same 
simulation code but with a softening length fixed in phys- 
ical units. The size evolution in this case is very similar 
to the simulations that uses a fixed comoving softening 
length. We also included various observational results 
which find slightly larger sizes at a given redshift but with 
a very similar evolution in time. On average there is a 
strong evolution in galaxy sizes: for the R\/2 oc (1 + z) a 
power law fit to all (dashed line) and only the quies- 
cent (solid line) simulated galaxies we find a value of 
a = -1.12 ± 0.13 and a = -1.44 ± 0.16, respectively. 
This is in good agreement with observed values - which 
are possibly a bit offset to higher values - and despite 
our statistical limitations we consider this trend robust. 
Possible simple explanations for an offset i n size within a 
semi-a nalytical framework are discussed in Shankar et al. 
(|20T1 . 

To demonstrate the physical origin for the size growth 
in the simulated galaxies we show in Fig. [2] the pro- 
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Fig. 1. — Left: Projected stellar half-mass radii of the simulated galaxies versus stellar masses for redshifts z=0 (blue circles), z=l (red 
triangles), z=2 (green triangles) and z=4 (purple stars). Quiescent galaxies with low specific star formation rates (sSFR < 0.3/%) have 
solid symbols, other galaxies have open sy mbols. Since z as 2 all galaxies evolve rapidly in size. The dashed line indicates the observed 
size-mass relation for early-type galaxies of (Nipoti ct al. 2009) with the one-sigma scatter indicated by the dotted lines. The z=0 mass cut 
of M* > 6.3 X 10 10 Mq for the galaxy sample is indicated by the vertical dashed line. Right: Projected stellar half-mass radii of galaxies 
with stellar masses Mr > 6.3 x 1O 1O M (see arrow on the left plot) as a function of redshift. The black symbols indicate the mean sizes 
at a given redshift with the error bars showing the standard deviation. The star forming galaxies (open symbols) and mean values are 
offset by 0.1 in redshift for clarity. The black lines show the result of a power law fit for all (dotted line) and quiescent (solid line) systems, 
respectively in good agreement with observed relations. Obs ervational estimates from different authors are given by the solid star symbols 
where the dott ed error bars show the ob served scatter (see [van Dokkum ct al. 2008; Ccnarro & Truiillo 2009) with the exeption of the 
observation by Ivan de Sande eTaT] (2011). Since this is a single object, here the error bars indicate the uncertainty of the measurement.. 
By z=3 all progenitor galaxies drop below our mass limit. 



jected half-mass radii presented in Fig. [T]at different red- 
shifts as a function of the ratio of stars accreted onto the 
galaxy Af^acc to the stars formed in-situ in the galaxy, 
M*,ins at the same redshifts. We consider a star parti- 
cle in the simulation as formed in-situ in the galaxy if 
it is created inside i? ga i(= 0.1 x i?™). Black bordered 
symbols indicate galaxies with stellar masses larger than 
6.3 x 1O 1O M0 whose size evolution is plotted in the right 
panel of Fig. Q] There is a clear correlation between 
the relative amount of ac creted stars and in-situ stars 
not only at redshift zero (jOser et al.ll2010t ) but also at 
high redshifts (z k, 2). This indicates that stellar ac- 
cretion drives the size evolution of the systems as soon 
as the accreted stars start to dominate the total mass 
(M*,acc/Mtsins > 1) at z » 2 as also predicted from semi- 
analytical modeling (jKhochfar fc Silk 2006b) . At earlier 
times the stell ar mass growth is dominated by in-situ 
star formation (|Oser et al.ll2010t ). i.e. the stars form out 
of cold gas that was able to radiate away a large frac- 
tion of its gravitational energy and thus leading to com- 
pact systems. The binding energy of the accreted stars, 
however, is retained and will increase the total energy 
content of the accreting galaxy, both by shock-heating 
the gas - which then can cool radiatively - as well as ex- 
panding the existing dar k matter and stellar components 
(|Johansson et a l. 2009b). This in general leads to more 
extended systems. 

In Fig. [3] we compare the density profiles of a subsam- 
ple of massive galaxies (M* > 2 x 1O 11 M ) at redshift 



2 and the present day. In agreement with iNaab et al.l 
(2009) we find that within the half-mass radius the high 
redshift systems are dominated by stars that formed in- 
situ while the contribution of accreted stars to the inner 
mass profile is small. At the present day the stellar mass 
inside the effective radius is dominated by accreted stars 
added at radii larger than > 1 kpc . This accretion is 
responsible for the strong size increase (lOser et alJl2010l) 
and is in agreement with the results from stacked imaging 
for massive galaxies at a constant number density that 
also show an increas e in surface densities pred ominantly 
in the outer regions |yan Dokk um et al.ll2010f ). 

4. REDSHIFT EVOLUTION OF VELOCITY DISPERSIONS 

In the left panel of Fig. [4] we show the central stellar 
line-of-sight velocity dispersions for the simulated galax- 
ies as a function of the stellar mass at redshift z=4 (pur- 
ple stars), z=2 (green triangles), z=l (red triangles) and 
z=0 (blue circles). The line-of-sight velocity dispersions 
have been calculated within 0.5 x R\/2 along the three 
principal axe s and then average d. The mass-dispersion 
relation from lNipoti et all (|2009f ) for galaxies more mas- 
sive than M* ~ 10 11 M Q is indicated by the dashed 
line. It is evident that at a given mass range the ve- 
locity dispersions of the galaxies systematically increase 
with redshift. This is illustrated in the right panel of 
Fig. 0] where we show the evolution of the projected ve- 
locity dispersion for galaxies with masses in the range 
of 6.3 x 10 n M o < M* < 1.6 x 10 12 Af Q (indicated by 
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Fig. 2. — Projected stellar half-mass radius of the simulated 
galaxies at different redshifts (see Fig. [TJ versus the fraction of 
stellar mass accreted (in major mergers, minor mergers and accre- 
tion events), M*, a cc, to the stellar mass formed in-situ, Af* i ns , in 
the galaxies. The black bordered symbols indicate systems more 
massive than M, > 6.3 X 1O 1O M0. At z > 2 galaxies with a higher 
fraction of accreted stars have larger sizes indicating that accretion 
of stellar systems drives the size evolution of massive galaxies. 




Fig. 3. — Stellar density as a function of radial distance to the 
galactic center at redshift (blue) and 2 (green) averaged over the 
most massive systems (M, > 2 X 10 11 Mq at z = 0). The arrows 
indicate the average half-mass radii at a given redshift. At z = 2 
we find that the stellar mass profile inside the half-mass radius is 
still dominated by stars that have formed in-situ (dotted lines). 
At z = the half-mass radii of our galaxies have significantly 
extended due to accreted stars (dashed lines) which dominate the 
total density (solid lines) at larger radii. 



the vertical lines in the left panel) as a function of red- 
shift. In this mass range the velocity dispersions drop 
from 262 ± 28kms _1 at z=2 to 177 ± 22kms _1 at z=0, 
a decrease of roughly a factor of 1.5. T he evolution is 
statist ically significant but weak (see also lHopkins et all 
2009bf). The black lines show a fit for the average ve- 
locity dispersions for all (dashed line) and the quies- 
cent (solid line) galaxies only. As for the sizes, we fit 
the redshift evolution of the velocity dispersions like 
01/2 oc (1 + z)p . Again we find a slightly stronger evo- 
lution for the quiescent systems (B = 0.44) than for all 
galaxi es in our samples (j3 — 0.37). Ivan de Sande et all 
(|2011h obtain a similar v alue of B = . 51 ±0. 07. Depend- 
ing on selection criteria iSaglia et al.l (|2010l ) find values 
for 8 ranging from 0.5 9 ± 0.1 to 0.19 ± 0.10. Follow- 
ing iCenarro fe Truiiilol (|2009f ) we compare to observa- 
tions of local ellipticals and measurements at higher red- 
shift (Idi Serego Alighieri et al.l 120051 : Ivan der Wei et al~l 
20051 120081: ICe narro fc Trujilldl2009( ) . In general we find 
a good agreement with the observations. 

5. STELLAR MERGER HISTORIES 



Similar to iCole et al.1 pOOOh . iNaab et al.l (|200l and 
iBezanson et al.l (|2009f ) demonstrated, using the virial 
theorem, how the size growth of a galaxy after a dis- 
sipationless merger event can be predicte d. Assuming 
energy conservation and p arabolic orbits (|Bensonl 120051 : 
iKhochfar fc BurkertJl2006j ) the final gravitational radius 
of the system is given by 



RgA (1 + V*) ' 



(1) 



where R g ^ and i? gj f are the initial and final gravitational 
radii (before and after the merger event) which we as- 
sume to be propor tional to the spherical half mass radii 
(|Naab et all 12009ft . Here, r\ = M a /Mi is the fractional 
mass increase during the merger and e = (v a 2 )/(vj 2 ) is 
the ratio of the mean square speeds of the accreted satel- 
lites and the initial system. This estimate implies that 
the accretion of weakly bound stellar systems (e <C 1) re- 
sults in a particularly efficient size increase. To test the 
validity of this simple estimate in our full cosmological 
simulation we follow the stellar merger histories of our 
simulated galaxies. We identify every satellite merging 
with the most massive progenitor of the central galaxy 
using a FOF-finder with a minimum number of 20 stellar 
particles (« 1.2 x 10 8 M©). At z — 2 all galaxies in our 
sample are more massive than » 10 .Mq, i.e. we re- 
solve mergers at least down to a mass-ratio of w 1 : 100. 
For the satellites as well as the host we compute the mass 
and the total internal velocity dispersion and with this in- 
formation and Eqn. [T]we estimate the size evolution since 
z=2. In Fig. [5] we show the estimated size growth for all 
galaxies exceeding the z=2 mass limit of 6.3 x 10 10 Af© 
since redshift 2 (green triangles) and 1 (red triangles), 
respectively. We compare this estimated size growth to 
the actual size growth that we directly measure from 
our simulated galaxies. Considering the simplifications 
used in Eqn. [1] - homologous, one-component systems 
merging on zero energy orbits - the predicted and actual 
grow th agree notably well. This confirms earlier find- 
ings ()Naab fc O striker 2009) that the simple formula is 
a good predictor even in a full cosmological context. This 
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Fig. 4. — Central (within 0.5 R1/2) projected velocity dispersion as a function of stellar mass at z=0 (blue circles), z=l (red triangles), 
z=2 (green triangles) and z=4 (purple stars). The relation for local galaxies from Nipoti ct al. (2009') are shown by the dashed line with 
the dotted lines indicating the scatter of the observed galaxies. At a given mass the velocity dispersion decreases significantly from z=4 
to z=0. The mass limits used for the right plot are indicated by the vertical dashed lines. Right: Central projected velocity dispersion 
of the simulated galaxies with masses in the range of 6.3 X 1O 1O M < Ah < 1.6 X 10 n A/ Q at any given redshift as a function of 
redshift. Solid symbols represent star forming galaxies and empty symbols show quie scent systems (offset by 0.1 in redshift for clarity) . 
Observational estimates from different authors are given by the solid star symbols (see Ccnarro & Truiillo 2009; van dc Sandc ct al. 2011; 
Martinez-Manso ct al. 2011) with the observed scatter given by the dotted error bars, where available. The black lines show the result of a 
power law fit for all (dashed line) and the quiescent (solid line) galaxies, respectively. The simulations indicate a mild dispersion evolution 
from ss 262 kms -1 at z=2 to es 177 kms -1 at z=0, in agreement with observations. 



15 



▼ size growth 


since z=l / 


A size growth 


since z = 2 - 
A / 




A 

A 

A/ 


▲ 


A 


















5 


10 15 



estimated size growth 



Fig. 5. — The size growth predicted by equation[T]in combination 
with the stellar merger histories compared to the actual size growth 
in the simulations of the galaxies more massive than M* = 6.3 X 
lO lo A/0 at z=2. The green triangles indicate the evolution between 
z=2 and z=0 the red triangles the evolution between z=l and z=0. 
The simple virial estimate is a good predictor for the actual size 
evolution. 



approximation however can only be valid if the assembly 
history for z < 2 is not dominated by dissipational pro- 
cesses wh ich is the case for the massive systems presented 
here (see lOser et a l. 2010 for the relevant analysis). 

To better understand the dominant assembly mecha- 
nism of our simulated massive galaxies we have computed 
the average mass-weighted merger mass-ratio for every 
galaxy since z=2. In the left panel of Fig. [S]we show the 
average mass-weighted merger mass-ratio as function of 
present-day stellar galaxy mass (blue dots). The depen- 
dence on galaxy mass is weak. The average values in 
bins of 0.2 dex in mass are shown by the black diamonds 
with one sigma errors bars. Overall the average mass- 
weighted merger mass-ratio is ~ 0.20 ±0.10. This makes 
'minor mergers' with mass ratios of 1:5 the dominant 
assembly m ode, on average, for the massive simulated 
galaxies feeelHirschma nrTet al.ll20TTl for a representation 
of a typical merger tree). The tendency of this ratio to 
change with the mass of the host system cannot be deter- 
mined by our calculations with any statistical certainty 
(the slope of the fitted curve is 0.05 ±0.18). However, we 
anticipate that for very low mass galaxies major mergers 
would become more important. If the slope of the mass 
function for satellites were <i(mN)/d(lnM) ~ —7, then 
the expectation would be that the mass- weighted merger 
ratio would be (2 — 7)/ (3 — 7) ~ 0.44 if dynamical friction 
were not a dominant process and (3 — 7)/(4 — 7) ~ 0.64 
if it were dominant. Thus, for low mass parent galax- 
ies, we would anticipate that the typical merger would 
be relatively 'major' with the ratio of parent to satellite 
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Fig. 6. — Left: The average mass-weighted stellar merger mass-ratios (since z=2) as a function of present-day galaxy mass (blue dots). 
The black diamonds show the binned averages within 0.2 dex in stellar mass with the one sigma error bars. Trends with galaxy mass are 
statistically not significant. The mass growth is dominated by minor mergers with a mass ratio of pa 1:5. Right: The average number- 
weighted merger mass-ratio (for all stellar mergers since z=2) as a function of present-day galaxy mass. There is a weak trend for more 
massive galaxies to experience relatively more minor mergers. On average most stellar mergers have mass-ratio of 1 : 16. 



being ^1:2. Here we note that a significant number 
of the simulated galaxies (7 out of 40, ~ 18 per cent) 
do not experience any merger with a mass ratio larger 
than 1:4, e.g. they have experienced no major merger 
since z=2 at all. Estimates of merger rates for mas- 
sive galaxies due to observations of disturbed system s 
(|Jogee et al.ll2009t iKavirai et aT1l2011tlLotz et al.ll201lT ). 
as wel l as semi-analyti c mode ls lead to similar results . 
E.g IKhochfar fc Silkl (I2009T) . IHopkins et all (|2010al) 
and lShankar et al.l ()2010al ) find, that massive early- type 
galaxies on average encounter less than one major dry 
merger since their formation epoch. This confirms pre- 
vious suggestions motivated by the dearth of compact 
galaxies in the nearby Universe, that a highly stochasti- 
cal process like major mergers can not be the main driver 
for the observed size evolution dBezanson "et~aT][200l 
iTruiillo et aT][2009l: iTaylor et al.ll2010D . However, major 
mergers do happen and will have an impact on the early- 
type galaxy population. They can contribute signifi- 
cantly to the final stellar m ass with minor mergers sti ll 
dominating the size growth (|Shankar et al.ll2010bl 1201 1ft . 
The observed merger rates, which are difficult to deter- 
mine, a re in the range of only ~ 1 major merger since 
z = 2 (|Bell et al.l [20061: iBluck et al.l [20091 INaab et al.l 
2006) which is consistent with our interpretation. 

In the right panel of Fig. [6] we show, the more conven- 
tionally defined, average number-weighted merger mass- 
ratio. The merger history since z=2 is clearly dominated 
by minor mergers with mass-ratios smaller than 1:10. 
Those mergers, however, do on average not add most of 
the mass to the systems. There is a slight trend for more 
massive galaxies to experience a larger relative number 
of minor mergers. Over the full mass range the average 
number-weighted merger mass ratio is ~ 0.062 ± 0.043, 



indicating that the typical merger was indeed very minor 
(1:16). 

6. CONCLUSION & DISCUSSION 

In this paper we use a sample of 40 cosmological re- 
simulations of individual massive galaxies to investigate 
the evolution of galaxy sizes and velocity dispersions with 
reds hift. The simulate d galaxies form in a two phase pro- 
cess (|Oser et al.ll2010D where the first phase at redshifts 
of z > 2 is dominated by a dissipative assembly. This 
formation phase is driven by in situ star formation re- 
sulting in compact galaxies having small sizes of r < 1.3 
kpc. The subsequent evolution of the galaxies at redshifts 
of z < 2 is do minated by accretion of stars in satellite 
stellar systems. INaab fe Ostrikerl (|2009f) and lOser et all 
( 2010) have shown that the accreted stellar systems pref- 
erentially settle into the outer parts of the galaxies, re- 
sulting in a gradual increase in their sizes until the sim- 
ulated galaxies closely follow the present-day mass-size 
relation. Between redshift 2 and our simulated galax- 
ies grow on average by a factor of ~ 5 — 6, whereas re- 
cent semi-ana lytical models find a smaller size increase 
of ~ 2-4 (IKhochfar fc Silkl 12006a!: IGuo et al.l [20TTI : 
iCovington et al.l 1201 If ) . At the present day 25 out of 
the 40 simulated galaxies are quiescent (sSFR < 0.3/in) 
and have structural parameters in agreement with ob- 
served local early-type galaxies. The underlying physical 
reason for the size growth for our simulated galaxies is 
stellar accretion, as can be seen in the strong positive 
correlation between the projected stellar half-mass radii 
and the fraction of accreted stellar material (Fig [2]). Our 
detailed analysis presented in this paper confirms that 
the stellar mate rial is predomina ntly accreted through 
minor mergers (jNaab et al.l 120091) , with typical galaxy 
mass-ratios of w 1 : 5. By number the merger history 
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is dominated by even more minor mergers with mass- 
ratios of » 1 : 16. A significant fraction (18 per cent) 
of the galaxies experience no major merger with mass- 
ratios larger than 1:4 since z=2 confirming previous sug- 
gestions, motivated by the lack of compact galaxies in the 
nearby Universe, that a highly stochastical process such 
as major mergers ca nnot be the main driver for the ob- 
served size evolution dBezanson et al.l2009tlTruiillo et al.l 
120091: iTavlor et al.l 12010ft . Semi-analytical models also 
find significant stellar mass growth due to minor mergers. 
These models, however, predict that for the most mas- 
sive galaxies m ajor mergers are becoming increasingly 
important ("e.g. iBaugh et alj|1996t iDe Lucia et"aLl l2006t 



IDe Lucia fc Blaizotl 120071: IGuo fc White! 12008ft . This is 
a result of the sharp drop-off in the galaxy mass func- 
tion due to AGN feedback, which is not followed in our 
simulations. 

For galaxies with masses above 6.3 x 1O 1O M0 our sim- 
ulated size evolution is in very good agreement (Fig [TJ 
with the observed size evolution of galaxi es with sim- 
ilar m asses at redshifts of z < 2 (e.g. iFranx et al.l 
2008). The evolution of the sizes can be well described 
by Ri/2 oc (1 + z) a with a = —1.12 for all galaxies and 
a = —1.44 for quiescent galaxies only. The size growth 
measured from the simulations is in good agreement with 
simple estimates from the virial theorem assuming en- 
ergy conservation during dissipationless merger events 
(|Naab et all 1 2559ft . 

The projected velocity dispersions for simulated galax- 
ies with masses around « 10 11 M Q decrease systemati- 
cally towards lower redshifts from s» 262 kms" 1 at z=2 
to 177 kms -1 at z = 0, again in good agr eement with 
observations (e.g. iCenarro fc Truiillol l2009h ■ Assuming 
an evolution as <7i/ 2 oc (1 + zf we find f3 = 0.37 for 
all galaxies and (3 = 0.44 for quiescent galaxies. Future 
observations might confirm this prediction. 

We conclude that in the absence of dissipation and as- 
sociated star formation a growth scenario dominated by 
minor stellar mergers, with less bound stars, is a viable 
physical process for explaining both the observed growth 
in size and the decrease in velocity dispersion of mas- 
sive early-type galaxies from z ~ 2 to the present-day. 
Accretion of systems not gravitationally bound to the 
central galaxy causes, as noted, substantial size growth. 
But it has another, dramatic, concomitant effect. As this 
mass becomes gravitationally bound, it releases a large 
amount of gravitational energy. T his process, which has 
been measured in our simula tions ([Khochfar fe Ostrikerl 
120081 : [Johansson et al.ll2009bft , and termed 'gravitational 
heating' can add ~ 10 59 5 ergs (i.e. ~ 10 43 erg/s) to the 
parent systems, causing heating of the ambient gas and 
reducing the central dark matter component. 

Despite these successes some obvious caveats concern- 



ing our simulations remain. Most importantly our simu- 
lated galaxies are overly efficie nt in transforming gas into 
stars ( Hirs chmann et af1l2011l ) and consequently the con- 
version efficiency of baryons into stars, even at z=2, in 
the massive galaxies in our simulated sample is overes- 
timated by roughly a factor of » 2 co mpared to pre- 
dictio ns from halo occupation models (Be hroozi et all 
120101 and references therein). This discrepancy is most 
probably due to the fact that our simulations neither 
include strong supernova-driven winds nor AGN feed- 
back from supermassive black holes and would be en- 
hanced if metal-line cooling was included. Observations 
and modeling have shown that strong galactic winds gen- 
erating significant outfl ows are ubiquitous at high red- 
shifts of z ~ 2 - 3 (e.g. iSteidel et all 120101: iGenel et all 
[2010ft and in our simulations this aspect is missing by 
construction. The effect of supernova driven winds and 
AGN feedback is differential with respect to the masses 
of galaxi es, with the former prima rily affecting smaller 
galaxies (jQppenheimer et ail 12010ft and the latter be- 
ing increasingly impor tant for more massive galaxies 
(jKormendv et al.l [20091 The proper inclusion of all the 
above mentioned physical effects would certainly lower 
the overall total stellar masses (both the in-situ and 
the accreted component). Still, the relatively simple 
two-phase formation scenario provides a viable model 
to physically explain the observed growth in size and 
decrease in velocity dispersion. However, this is an is- 
sue clearly deserv ing further studies on the effect of 
AGN f eedback (e.glSazonov et alJl200llDi Matteo et all 
20051: iSpringelet all 120051: I Johansson et al.n2 009a; 



Booth fc Schavel2009tlDebuhr et al. 201ft lOstriker' et all 



2010t iSchave et aill2010ft. r a diativ e feedback from stars 
(e.g [Qppenheimer fc Pavel 120061 : iPetkova fc Springel 
20091 : iHopkins et alll201lT) and feedback from supernovae 
type II (e.g iScannapieco et alJ 12008: Sch ave et al.l 1 20101 
Sales et all 12010ft and la Ce.g iCi otti fc OstrikeiT 120071: 
Qppenheimer fc Pavel 120081 Preferentially this will be 
investigated with the help of a large sample of zoom sim- 
ulations with better statistics, as presented here, with 
the aim of studying how these processes would affect in 
detail the resulting size growth and velocity dispersion 
evolution of massive galaxies. 



We thank Ignacio Trujillo, Pieter van Dokkum, Simon 
P. M. White, and Shardha Jogee as well as the anony- 
mous referee for very useful comments on the manuscript. 
Part of the simulations were performed at the Princeton 
PICSciE HPC center. This research was supported by 
the PFG cluster of excellence 'Origin and Structure of 
the Universe' as well as the PFG priority program 1177. 



REFERENCES 



Arnold, J. A., Romanowsky, A. J., Brodie, J. P., Chomiuk, L., 
Spitler, L. R., Strader, J., Benson, A. J., & Forbes, D. A. 2011, 
ApJ, 736, L26+ 

Baugh, C. M., Cole, S., & Frenk, C. S. 1996, MNRAS, 283, 1361 
Behroozi, P. S., Conroy, C, & Wechsler, R. H. 2010, ApJ, 717, 
379 



Bell, E. F., Naab, T., Mcintosh, D. H., Somerville, R. S., 

Caldwell, J. A. R., Barden, M., Wolf, C, Rix, H.-W., Beckwith, 
S. V., Borch, A., Haussler, B., Heymans, C, Jahnke, K., Jogee, 
S., Koposov, S., Meisenheimer, K., Peng, C. Y., Sanchez, S. F., 
& Wisotzki, L. 2006, ApJ, 640, 241 

Benson, A. J. 2005, MNRAS, 358, 551 

Bczanson, R., van Dokkum, P. C, Tal, T., Marchesini, D., Kriek, 
M., Franx, M., & Coppi, P. 2009, ApJ, 697, 1290 



Size and dispersion evolution of early-type galaxies 



9 



Bluck, A. F. L., Conselice, C. J., Bouwens, R. J., Daddi, E., 
Dickinson, M., Papovich, C, & Yan, H. 2009, MNRAS, 394, 
L51 

Booth, C. M. & Schaye, J. 2009, MNRAS, 398, 53 

Bournaud, F., Chapon, D., Teyssier, R., Powell, L. C, Elmcgrccn, 
B. G., Elmegreen, D. M., Due, P.-A., Contini, T., Epinat, B., & 
Shapiro, K. L. 2011, ApJ, 730, 4 

Buitrago, F., Trujillo, I., Conselice, C. J., Bouwens, R. J., 
Dickinson, M., & Yan, H. 2008, ApJ, 687, L61 

Cappellari, M., di Scrcgo Alighicri, S., Cimatti, A., Daddi, E., 
Renzini, A., Kurk, J. D., Cassata, P., Dickinson, M., 
Franccschini, A., Mignoli, M., Pozzetti, L., Rodighiero, G., 
Rosati, P., & Zamorani, G. 2009, ApJ, 704, L34 

Carrasco, E. R., Conselice, C. J., & Trujillo, I. 2010, MNRAS, 

405, 2253 

Cassata, P., Giavalisco, M., Guo, Y., Renzini, A., Ferguson, H., 
Koekemoer, A. M., Salimbeni, S., Scarlata, C, Grogin, N. A., 
Conselice, C. J., Dahlen, T., Lotz, J. M., Dickinson, M., & Lin, 
L. 2011, ArXiv e-prints, 1106.4308 

Cenarro, A. J. & Trujillo, I. 2009, ApJ, 696, L43 

Cimatti, A., Cassata, P., Pozzetti, L., Kurk, J., Mignoli, M., 
Renzini, A., Daddi, E., Bolzonella, M., Brusa, M., Rodighiero, 
G., Dickinson, M., Franccschini, A., Zamorani, G., Berta, S., 
Rosati, P., & Halliday, C. 2008, A&A, 482, 21 

Ciotti, L. & Ostriker, J. P. 2007, ApJ, 665, 1038 

Coccato, L., Gerhard, O., Arnaboldi, M., & Ventimiglia, G. 2011, 
ArXiv e-prints, 1108.3834 

Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, 
MNRAS, 319, 168 

Covington, M. D., Primack, J. R., Porter, L. A., Croton, D. J., 
Somerville, R. S., & Dckel, A. 2011, ArXiv c-prints 

Daddi, E., Renzini, A., Pirzkal, N., Cimatti, A., Malhotra, S., 
Stiavelli, M., Xu, C, Pasquali, A., Rhoads, J. E., Brusa, M., di 
Serego Alighicri, S., Ferguson, H. C, Koekemoer, A. M., 
Moustakas, L. A., Panagia, N., & Windhorst, R. A. 2005, ApJ, 
626, 680 

Damjanov, I., Abraham, R. G., Glazcbrook, K., McCarthy, P. J., 
Caris, E., Carlbcrg, R. G., Chen, H.-W., Crampton, D., Green, 

A. W., j0rgensen, I., Juneau, S., Le Borgne, D., Marzke, R. O., 
Mcntuch, E., Murowinski, R., Roth, K., Savaglio, S., & Yan, H. 
2011, ArXiv c-prints, 1108.0656 

Damjanov, I., McCarthy, P. J., Abraham, R. G., Glazebrook, K., 
Yan, H., Mentuch, E., Le Borgne, D., Savaglio, S., Crampton, 
D., Murowinski, R., Juneau, S., Carlbcrg, R. G., j0rgensen, I., 
Roth, K., Chen, H., & Marzke, R. O. 2009, ApJ, 695, 101 

Dc Lucia, G. & Blaizot, J. 2007, MNRAS, 375, 2 

De Lucia, G., Springel, V., White, S. D. M., Croton, D., & 
Kauffmann, G. 2006, MNRAS, 366, 499 

Debuhr, J., Quataert, E., Ma, C, & Hopkins, P. 2010, MNRAS, 

406, L55 

Dekel, A., Birnboim, Y., Engel, G., Freundlich, J., Goerdt, T., 
Mumcuoglu, M., Neistein, E., Pichon, C, Teyssier, R., & 
Zinger, E. 2009, Nature, 457, 451 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 
604 

di Scrcgo Alighieri, S., Vernet, J., Cimatti, A., Lanzoni, B., 
Cassata, P., Ciotti, L., Daddi, E., Mignoli, M., Pignatelli, E., 
Pozzetti, L., Renzini, A., Rettura, A., & Zamorani, G. 2005, 
A&A, 442, 125 

Faber, S. M., Willmer, C. N. A., Wolf, C, Koo, D. C, Weiner, 

B. J., Newman, J. A., Im, M., Coil, A. L., Conroy, C, Cooper, 
M. C, Davis, M., Finkbeiner, D. P., Gerke, B. F., Gebhardt, 
K, Groth, E. J., Guhathakurta, P., Harker, J., Kaiser, N., 
Kassin, S., Klcinhcinrich, M., Konidaris, N. P., Kron, R. G., 
Lin, L., Luppino, G., Madgwick, D. S., Meisenhcimcr, K., 
Noeske, K. G., Phillips, A. C, Sarajcdini, V. L., Schiavon, 

R. P., Simard, L., Szalay, A. S., Vogt, N. P., & Yan, R. 2007, 
ApJ, 665, 265 

Fan, L., Lapi, A., Bressan, A., Bernardi, M., De Zotti, G., & 

Danese, L. 2010, ApJ, 718, 1460 
Fan, L., Lapi, A., Dc Zotti, G., & Danese, L. 2008, ApJ, 689, L101 
Fcldmann, R., Carollo, C. M., Mayer, L., Renzini, A., Lake, G., 

Quinn, T., Stinson, G. S., & Yepes, G. 2010, ApJ, 709, 218 
Forbes, D. A., Spitler, L. R., Strader, J., Romanowsky, A. J., 

Brodic, J. P., & Foster, C. 2011, MNRAS, 413, 2943 
Franx, M., van Dokkum, P. G., Schreiber, N. M. F., Wuyts, S., 

Labbe, I., & Toft, S. 2008, ApJ, 688, 770 



Genel, S., Naab, T., Gcnzel, R., Forster Schreiber, N. M., 
Sternberg, A., Oser, L., Johansson, P. H., Dave, R., 
Oppcnhcimer, B. D., & Burkcrt, A. 2010, ArXiv c-prints 

Guo, Q., White, S., Boylan-Kolchin, M., De Lucia, G., 

Kauffmann, G., Lemson, G., Li, C, Springel, V., & Weinmann, 
S. 2011, MNRAS, 413, 101 

Guo, Q., White, S., Li, C, & Boylan-Kolchin, M. 2010, MNRAS, 
404, 1111 

Guo, Q. & White, S. D. M. 2008, MNRAS, 384, 2 

Guo, Y., Mcintosh, D. H., Mo, H. J., Katz, N., van den Bosch, 

F. C, Weinberg, M., Weinmann, S. M., Pasquali, A., & Yang, 

X. 2009, MNRAS, 398, 1129 
Haardt, F. & Madau, P. 1996, ApJ, 461, 20 

Hirschmann, M., Naab, T., Somerville, R., Burkcrt, A., & Oser, 

L. 2011, ArXiv e-prints, 1104.1626 
Hopkins, P. F., Bundy, K., Croton, D., Hernquist, L., Keres, D., 

Khochfar, S., Stewart, K., Wetzel, A., & Younger, J. D. 2010a, 

ApJ, 715, 202 

Hopkins, P. F., Bundy, K., Hernquist, L., Wuyts, S., & Cox, T. J. 

2010b, MNRAS, 401, 1099 
Hopkins, P. F., Bundy, K., Murray, N., Quataert, E., Lauer, 

T. R., & Ma, C.-P. 2009a, MNRAS, 398, 898 
Hopkins, P. F., Hernquist, L., Cox, T. J., Keres, D., & Wuyts, S. 

2009b, ApJ, 691, 1424 
Hopkins, P. F., Quataert, E., & Murray, N. 2011, ArXiv e-prints 
Hyde, J. B. & Bernardi, M. 2009, MNRAS, 394, 1978 
Jogee, S., Miller, S. H., Penner, K., Skelton, R. E., Conselice, 
C. J., Somerville, R. S., Bell, E. F., Zheng, X. Z., Rix, H.-W., 
Robaina, A. R., Barazza, F. D., Barden, M., Borch, A., 
Bcckwith, S. V. W., Caldwell, J. A. R., Peng, C. Y., Heymans, 
C, Mcintosh, D. H., HauBler, B., Jahnke, K., Meisenheimer, 
K., Sanchez, S. F., Wisotzki, L., Wolf, C, & Papovich, C. 2009, 
ApJ, 697, 1971 

Johansson, P. H., Naab, T., & Burkert, A. 2009a, ApJ, 690, 802 
Johansson, P. H., Naab, T., & Ostriker, J. P. 2009b, ApJ, 697, 
L38 

Joung, M. R., Cen, R., & Bryan, G. L. 2009, ApJ, 692, LI 
Kauffmann, G. 1996, MNRAS, 281, 487 

Kaviraj, S., Tan, K.-M., Ellis, R. S., & Silk, J. 2011, MNRAS, 
411, 2148 

Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, MNRAS, 
363, 2 

Khochfar, S. & Burkert, A. 2006, A&A, 445, 403 
Khochfar, S. & Ostriker, J. P. 2008, ApJ, 680, 54 
Khochfar, S. & Silk, J. 2006a, ApJ, 648, L21 
— . 2006b, MNRAS, 370, 902 
— . 2009, MNRAS, 397, 506 

Kormendy, J., Fisher, D. B., Cornell, M. E., & Bender, R. 2009, 
ApJS, 182, 216 

Kriek, M., van der Wei, A., van Dokkum, P. G., Franx, M., & 

Illingworth, G. D. 2008, ApJ, 682, 896 
Kriek, M., van Dokkum, P. G., Franx, M., Quadri, R., Gawiser, 

E. , Herrera, D., Illingworth, G. D., Labbe, I., Lira, P., 
Marchesini, D., Rix, H.-W., Rudnick, G., Taylor, E. N., Toft, 
S., Urry, C. M., & Wuyts, S. 2006, ApJ, 649, L71 

Longhetti, M., Saracco, P., Scvcrgnini, P., Delia Ceca, R., 

Mannucci, F., Bender, R., Drory, N., Feulner, G., & Hopp, U. 

2007, MNRAS, 374, 614 
Lotz, J. M., Jonsson, P., Cox, T. J., Croton, D., Primack, J. R., 

Somerville, R. S., & Stewart, K. 2011, ArXiv e-prints, 

1108.2508 

Martincz-Manso, J., Guzman, R., Barro, G., Cenarro, J., 

Perez-Gonzalez, P., Sanchcz-Blazquez, P., Trujillo, I., Balcells, 
M., Cardiel, N., Gallego, J., Hempel, A., & Prieto, M. 2011, 
ApJ, 738, L22+ 

Moslch, M., Williams, R. J., Franx, M., & Kriek, M. 2011, ApJ, 
727, 5 

Moster, B. P., Somerville, R. S., Maulbetsch, C, van den Bosch, 

F. C, Maccio, A. V., Naab, T., & Oser, L. 2010, ApJ, 710, 903 
Naab, T., Johansson, P. H., & Ostriker, J. P. 2009, ApJ, 699, 

L178 

Naab, T., Johansson, P. H., Ostriker, J. P., & Efstathiou, G. 

2007, ApJ, 658, 710 
Naab, T., Khochfar, S., & Burkert, A. 2006, ApJ, 636, L81 
Naab, T. & Ostriker, J. P. 2009, ApJ, 690, 1452 
Newman, A. B., Ellis, R. S., Trcu, T., & Bundy, K. 2010, ApJ, 

717, L103 



10 



Oser et al. 



Nipoti, C, Londrillo, P., & Ciotti, L. 2003, MNRAS, 342, 501 
Nipoti, C, Treu, T., Auger, M. W., & Bolton, A. S. 2009, ApJ, 
706, L86 

Oppenheimer, B. D. & Dave, R. 2006, MNRAS, 373, 1265 
— . 2008, MNRAS, 387, 577 

Oppenheimer, B. D., Dave, R., Keres, D., Fardal, M., Katz, N., 

Kollmeier, J. A., & Weinberg, D. H. 2010, MNRAS, 406, 2325 
Oser, L., Ostriker, J. P., Naab, T., Johansson, P. H., & Burkcrt, 

A. 2010, ApJ, 725, 2312 
Ostriker, J. P., Choi, E., Ciotti, L., Novak, G. S., & Proga, D. 

2010, ApJ, 722, 642 
Petkova, M. & Springel, V. 2009, MNRAS, 396, 1383 
Power, O, Navarro, J. F., Jenkins, A., Frenk, C. S., White, 

S. D. M., Springel, V., Stadel, J., & Quinn, T. 2003, MNRAS, 

338, 14 

Ricciardelli, E., Trujillo, I., Buitrago, F., & Conselice, C. J. 2010, 
MNRAS, 406, 230 

Ryan, Jr., R. E., McCarthy, P. J., Cohen, S. H., Yan, H., Hathi, 
N. P., Koekemoer, A. M., Rutkowski, M. J., Mechtley, M. R., 
Windhorst, R. A., O'Connell, R. W., Balick, B., Bond, H. E., 
Bushouse, H., Calzetti, D., Crockett, R. M., Disney, M., 
Dopita, M. A., Frogel, J. A., Hall, D. N. B., Holtzman, J. A., 
Kaviraj, S., Kimble, R. A., MacKenty, J., Mutchler, M., 
Paresce, F., Saha, A., Silk, J. I., Trauger, J., Walker, A. R., 
Whitmore, B. O, & Young, E. 2010, ArXiv e-prints 

Saglia, R. P., Sanchcz-Blazquez, P., Bender, R., Simard, L., 

Desai, V., Aragon-Salamanca, A., Milvang- Jensen, B., Halliday, 

0, Jablonka, P., Noll, S., Poggianti, B., Clowe, D. I., De Lucia, 
C, Pello, R., Rudnick, C, Valentinuzzi, T., White, S. D. M., & 
Zaritsky, D. 2010, A&A, 524, A6+ 

Sales, L. V., Navarro, J. F., Schaye, J., Vecchia, C. D., Springel, 

V., & Booth, C. M. 2010, MNRAS, 409, 1541 
Saracco, P., Longhctti, M., & Andreon, S. 2009, MNRAS, 392, 

718 

Sazonov, S. Y., Ostriker, J. P., & Sunyaev, R. A. 2004, MNRAS, 
347, 144 

Scannapieco, O, Tissera, P. B., White, S. D. M., k, Springel, V. 

2008, MNRAS, 389, 1137 
Schaye, J., Dalla Vecchia, O, Booth, C. M., Wiersma, R. P. O, 

Theuns, T., Haas, M. R., Bertone, S., Duffy, A. R., McCarthy, 

1. C, & van de Voort, F. 2010, MNRAS, 402, 1536 
Shankar, F., Marulli, F., Bernardi, M., Boylan-Kolchin, M., Dai, 

X., & Khochfar, S. 2010a, MNRAS, 405, 948 

Shankar, F., Marulli, F., Bernardi, M., Dai, X., Hyde, J. B., & 
Sheth, R. K. 2010b, MNRAS, 403, 117 

Shankar, F., Marulli, F., Bernardi, M., Mei, S., Meert, A., & 
Vikram, V. 2011, ArXiv e-prints, 1105.6043 

Spergel, D. N., Bean, R., Dore, O., Nolta, M. R., Bennett, C. L., 
Dunkley, J., Hinshaw, G., Jarosik, N., Komatsu, E., Page, L., 
Peiris, H. V., Verde, L., Halpern, M., Hill, R. S., Kogut, A., 
Limon, M., Meyer, S. S., Odcgard, N., Tucker, G. S., Weiland, 
J. L., Wollack, E., & Wright, E. L. 2007, ApJS, 170, 377 

Springel, V. 2005, MNRAS, 364, 1105 

Springel, V., Di Matteo, T., & Hernquist, L. 2005, ApJ, 620, L79 

Springel, V. & Hernquist, L. 2003, MNRAS, 339, 289 

Steidel, C. C, Erb, D. K., Shapley, A. E., Pettini, M., Reddy, N., 

Bogosavljevic, M., Rudie, G. C, & Rakic, O. 2010, ApJ, 717, 

289 

Szomoru, D., Franx, M., van Dokkum, P. G., Trenti, M., 

Illingworth, G. D., Labbe, I., Bouwens, R. J., Oesch, P. A., & 
Carollo, C. M. 2010, ApJ, 714, L244 



Taylor, E. N., Franx, M., Glazebrook, K., Brinchmann, J., van 
dcr Wei, A., & van Dokkum, P. G. 2010, ApJ, 720, 723 

Toft, S., van Dokkum, P., Franx, M., Labbe, I., Forster Schreiber, 
N. M., Wuyts, S., Webb, T., Rudnick, G., Zirm, A., Kriek, M., 
van der Werf, P., Blakeslee, J. P., Illingworth, G., Rix, H.-W., 
Papovich, C, & Moorwood, A. 2007, ApJ, 671, 285 

Trujillo, I., Ccnarro, A. J., de Lorenzo-Caceres, A., Vazdekis, A., 
de la Rosa, I. G., & Cava, A. 2009, ApJ, 692, L118 

Trujillo, I., Conselice, C. J., Bundy, K., Cooper, M. C, 
Eisenhardt, P., & Ellis, R. S. 2007, MNRAS, 382, 109 

Trujillo, I., Ferreras, I., & de La Rosa, I. G. 2011, MNRAS, 415, 
3903 

Trujillo, I., Forster Schreiber, N. M., Rudnick, G., Barden, M., 

Franx, M., Rix, H.-W., Caldwell, J. A. R., Mcintosh, D. H., 

Toft, S., Haussler, B., Zirm, A., van Dokkum, P. G., & Labbe, 

I. 2006, ApJ, 650, 18 
van de Sande, J., Kriek, M., Franx, M., van Dokkum, P. G., 

Bczanson, R., Whitaker, K. E., Brammer, G., Labbe, I., Groot, 

P. J., & Kaper, L. 2011, ApJ, 736, L9+ 
van der Wei, A., Bell, E. F., van den Bosch, F. C, Gallazzi, A., & 

Rix, H. 2009, ApJ, 698, 1232 
van der Wei, A., Franx, M., van Dokkum, P. G., Rix, H.-W., 

Illingworth, G. D., & Rosati, P. 2005, ApJ, 631, 145 
van der Wei, A., Holden, B. P., Zirm, A. W., Franx, M., Rettura, 

A., Illingworth, G. D., & Ford, H. C. 2008, ApJ, 688, 48 
van der Wei, A., Rix, H.-W., Wuyts, S., McGrath, E. J., 

Koekemoer, A. M., Bell, E. F., Holden, B. P., Robaina, A. R., 

& Mcintosh, D. H. 2011, ApJ, 730, 38 
van Dokkum, P. G. & Brammer, G. 2010, ApJ, 718, L73 
van Dokkum, P. G., Franx, M., Kriek, M., Holden, B., Illingworth, 

G. D., Magee, D., Bouwens, R., Marchesini, D., Quadri, R., 

Rudnick, G., Taylor, E. N., & Toft, S. 2008, ApJ, 677, L5 
van Dokkum, P. G., Kriek, M., & Franx, M. 2009, Nature, 460, 

717 

van Dokkum, P. G., Quadri, R., Marchesini, D., Rudnick, G., 
Franx, M., Gawiser, E., Herrera, D., Wuyts, S., Lira, P., Labbe, 
I., Maza, J., Illingworth, G. D., Forster Schreiber, N. M., 
Kriek, M., Rix, H., Taylor, E. N., Toft, S., Webb, T., & Yi, 
S. K. 2006, ApJ, 638, L59 

van Dokkum, P. G., Whitaker, K. E., Brammer, G., Franx, M., 
Kriek, M., Labbe, I., Marchesini, D., Quadri, R., Bezanson, R., 
Illingworth, G. D., Muzzin, A., Rudnick, G., Tal, T., & Wake, 
D. 2010, ApJ, 709, 1018 

Weinzirl, T., Jogee, S., Conselice, C. J., Papovich, C, Chary, 
R. R., Bluck, A. F. L., Gruetzbauch, R., Buitgrago, F., 
Mobasher, B., Lucas, R. A., Dickinson, M., & Bauer, A. E. 
2011, ArXiv e-prints, 1107.2591 

Williams, R. J., Quadri, R. F., Franx, M., van Dokkum, P., & 
Labbe, I. 2009, ApJ, 691, 1879 

Williams, R. J., Quadri, R. F., Franx, M., van Dokkum, P., Toft, 
S., Kriek, M., & Labbe, I. 2010, ApJ, 713, 738 

Wuyts, S., Cox, T. J., Hayward, C. C, Franx, M., Hernquist, L., 
Hopkins, P. F., Jonsson, P., & van Dokkum, P. G. 2010, ApJ, 
722, 1666 

Zirm, A. W., van der Wei, A., Franx, M., Labbe, I., Trujillo, I., 
van Dokkum, P., Toft, S., Daddi, E., Rudnick, G., Rix, H., 
Rottgering, H. J. A., & van der Werf, P. 2007, ApJ, 656, 66 



